Leading-order behavior of the correlation energy in the uniform electron gas 
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We show that, in the high-density limit, restricted M0ller-Plesset (RMP) perturbation theory 
yields -Ermp = 7r~ 2 (l — In 2) lnr a + O (r°) for the correlation energy per electron in the uniform 
electron gas, where r s is the Seitz radius. This contradicts an earlier derivation which yielded 
-Ermp = 0(ln |lnr s |). The reason for the discrepancy is explained. 
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We consider a paramagnetic system of N interacting 
electrons confined in a cubic box, with edges of length L 
and volume f2 = L 3 . We also assume a uniform charge 
density background of density equal in magnitude, but 
opposite in sign, to the average electron density p = N/il. 
In the thermodynamic limit, both N and VL tend to infin- 
ity in such a way that the system becomes homogeneous 
with a uniform density p, related to the Seitz radius by 
the relation r s = (4np/3) , and is often called jellium 

an. 

It is convenient to consider a reduced Hamiltonian 
(i.e. one that is scaled by the number of electrons) and, 
in atomic units, this is 
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N- 1 T + H P . 



Hi 



b-b 



where the operator 
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corresponds to the kinetic energy of the electrons, and 

N 



Y — 

N 



dR 



-ffb-b = — 



n - R\ 
dR 1 dR 2 



(3) 

(4) 
(5) 



represent the electron-electron, electron-background and 
background-background interactions, respectively. [|| 
The term 7?b-b is a known constant Q and may be ig- 
nored. 

In perturbation theory^, we introduce a partition 



H = H {0) + V, 



(6) 
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where the perturbation V is assumed small (in some 
sense) compared to the zeroth-order Hamiltonian H^ '. 
This yields an expansion of the (reduced) energy 

E = eM+eM+E<V + .... (7) 
The zeroth-, first- and second-order energies are given by 



$(o,/) 



ff(0) 



(0,«) 



(8) 
(9) 

(10) 



where £ (0 < 0) = £ (0) and *(°<°) = ¥°\ and the zeroth- 
order ground state (£ = 0) and excited states (£ > 0) 
wave functions satisfy 
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There are many ways to partition H but not all are 
equally effective. In this Brief Report, we will consider 
three: the non-interacting (NI), restricted M0ller-Plesset 
H] (RMP), and unrestricted M0ller-Plesset (UMP) parti- 
tions. 

If we adopt the NI partition, we have 



MO) _ rf, 



Vni = H c . c + H c _b . 
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The zeroth-order wave functions are Slater deter- 



minants of plane- wave orbitals 



tpk (r) = — exp (i k ■ r) . 



(13) 



with orbital energies 



e N i(fc) 



The l-th excited determinant "f^i has the energy 
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where the sum over k takes into account all the plane 

waves used to build ^ni^> *- e - a ^ the occupied orbitals 
in the state £. For the special case I = 0, all the orbitals 
up to the Fermi level are occupied. 

Introducing a = (97r/4) 1 ^ 3 , one finds 0-Q that 
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which are the kinetic and exchange energies, respectively. 

Unfortunately, although the correlation energy [1JJ 

E C = E- £(°> - (17) 

of jellium is known fl2| to be finite for any r s > 0, the 
second-order energy Eq. ([TO]) is infinite. However, the 
leading-order contribution can be extracted from Eq. (|10[) 
and, henceforth, we will use E^ to refer to that contri- 
bution. 

After transforming into momentum space and scaling 
the momenta by the wave vector ftp = a/r s so that the 

Fermi sphere has unit radius, one finds [HI, [l4[ that E^ 
consists of a direct ( "ring-diagram" ) contribution 
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and an exchange contribution 
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(19) 



32tt 5 JJJ q 2 \q + k x -fc 2 | 2 Ae NI ' 

In these integrals, the excitation vector q has the domain 

p < \q\ < oo, (20) 

where (3 oc ^/rj [l5| . and the occupied-orbital vectors k\ 
and fc 2 have domains 



|fci+g|>l, (21) 
\k 2 -q\>l, (22) 



Ifcil < 1, 
|fc 2 | < 1, 



The lower bound for q in Eq. (|20|) is due to the screening 
effect of the Coulomb field by the collective electron mo- 
tions, and can be derived using the plasma theory of the 
free-electron gas (l5l - [l9| . The orbital energy difference is 

Ae N i = e NI (|fei + q\) + e N i (|fe 2 - q\) 

- em(fa) - e N i(fc 2 ). (23) 

The exchange term E^ h ^ is finite [2(J and, for small r s , 
is dominated by the ring-diagram term 



(24) 



which Macke showed [l3[ to depend logarithmically on 
r s . One may wonder, however, whether this logarithmic 



term arises when the Hamiltonian is partitioned differ- 
ently mi. 

If we adopt the RMP partition [22j], we have 
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RMP 
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= ^F(r 2 ), Vrmp = H-J2F{ri), (25) 
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where the Fock operator defined by 
1. 
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dr 2 (26) 



includes kinetic and exchange terms but not Hartrec 
terms because of their cancelation by the -ff e _b term. 

The RMP zeroth-order wave functions fjj^p are again 
determinants of plane- wave orbitals (|13[) . but the orbital 
energies are now different and it can be shown [23], [24| 
that 



CRMp(fc) = CNi(fc) -ex(fe)) 

cot 



The additional term 



1 - ft 2 
ex(fc) = l + ^^ln 



1 + ft 



1 - ft 



arises from the exchange terms in Eq. (|26j) . Thus, 



RMP — *NI ' 



but 
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The zeroth- and first-order energies are now given by 
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and comparing Eqs (JT6J) and ([31]) reveals the important 
relation 



E 
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NI 



^NI - A RMP + fc RMP - -&RHF, (32) 
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where -Erhf is the reduced RHF energy. 

(2) 

The ring-diagram contribution to E^J, p is 
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which differs from Eq. (| only by the denominator 



Acrmp — AeNi Aex, 

cot 



(34) 



where 



Ae x = ex (|fci + q\) + ex (|fc a - g|) 

- ex(fti) - e x (ft 2 ). (35) 
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The behavior of Ep^ p is dominated [24| by contributions 
in the neighborhood of the Fermi sphere (i.e. q « 0). On 
the domains (j2"Tj) and (|22j) . we have 
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Therefore, we have 



Ae RM p « u + v (u In - + v In - , 

an \ 2 2/ 



where we have introduced 
fci • q 



k 2 ■ q 



ki fc 2 
Substituting (|38|) into ([33| and using the relations 

dki = 2nk l sin g?# dA^ — du dki , 
|fei + q\ > 1 ^ 1 - u < ki < 1, 
(with similar expressions for k 2 ) then yields 
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Since the most important contribution comes from small q, we have set the upper bound of the integral (|42 
convenient value of unity. Expanding for small r s and integrating over q yields 
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(42) 



to a 



(43) 



which is identical, in the high-density (i.e. small-r s ) limit, 
to E^ a \ The present result can also be obtained from 
(|42j) by first switching to polar coordinates (u = rcosO 
and v = rsin#), integrating over the radial part, care- 
fully taking the r s — > limit, and finally performing the 
remaining angular integration. The latter derivation rig- 
orously justifies the small- r s expansion [25| . 

In a similar investigation more than 20 years ago [2lj , 
Handler claimed to show that 

4MP = 0(ln|lnr s |). (44) 

f2) 

This claim, which implies that Ep^p grows more slowly 

(2) 

with r s than E^i , obviously disagrees with our result in 
Eq. ([43]). However, in his analog of Eq. (g2]), Handler 
drops the u + v term and ignores the r s /an factor [2(3]. 
The fact that j3 oc ^fr~ s means that Handler's neglect of 
the u + v term is incorrect. 

(2) . (2) 

It may be surprising that Ep^p is the same as Ejff, 

because Hpyip seems a better starting point than . 
However, this is not the first time that the RHF treat- 
ment of jellium has been disappointing. For example, the 
RHF bandwidth, e(l) — e(0), is greater than the NI band- 
width, which disagrees with experiments on simple met- 
als, where a small reduction is observed [27], [28[. More- 
over, the logarithmic dependence of the eigenvalues (|2"T|) 



leads to a divergent derivative of CRMp(fc) at the surface 
of the Fermi sphere (k = 1) and this leads to incorrect 
dependence of the electronic specific heat on tempera- 
ture. Experimentally, a linear dependence with a prefac- 
tor close to the NI value is observed [29|, [3(| . 

One may hope that a different, and superior, perturba- 
tion series can be obtained by adopting the UMP parti- 
tion, that is, by using the UHF wavefunction of jellium as 
the star ting point. After all, as Overhauser showed long 
ago [3l|, l32j |. the RHF solution of jellium is unstable with 
respect to a lower-energy UHF solution, for all r s [33j . 
However, we expect that there will be serious issues with 
the convergence of the UMP perturbation series [H, [35j] 
and we have not considered this alternative in detail. 

In conclusion, we have shown that the correlation en- 

(2) 

ergy Ep^ p from RMP perturbation theory, i.e. using a 

(2) 

RHF starting point, is the same as the E^ from con- 
ventional NI perturbation theory. Although it is nearly 

impossible to test experimentally this result, this cor- 

(2) 

rects an earlier study which claimed that Ep^p is sub- 
logarithmic. 
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